arXiv: 1503.08170v2 [cond-mat.mes-hall] 17Jun2015 


Cooling of Levitated Graphene Nanoplatelets in High Vacuum 


Pavel Nagornykh and Joyce E. Coppock 
Department of Physics, University of Maryland, College Park, Maryland 20742, USA 

B. E. Kane 

Joint Quantum Institute, University of Maryland, College Park, Maryland 20742, USA and 
Laboratory for Physical Sciences, College Park, Maryland 20740, USA 

We demonstrate cooling of the center of mass motion of charged graphene nanoplatelets levitated 
in a quadrupole ion trap in high vacuum down to temperatures of 20 K. Parametric feedback based 
on optical measurements of particle motion was used to achieve the particle cooling at pressure 
p < 10 -6 Torr, and cooling along all three axes of motion was observed. Dependence of cooling 
on the electric fields was measured by varying DC voltages on a set of auxiliary electrodes used 
to spatially shift the trap minimum. Methods to calibrate mass and charge of the nanoplatelet by 
measuring its motion frequency dependence on discharge were also explored. 


While trapping and cooling technologies have been 
primarily applied to atomic systems m, they have also 
found applications in the study of nanoscale systems of 
condensed matter @hhi. One of many 2D materials 
that could benefit from these techniques is graphene, 
a material that garnered attention after its discovery 
in 2004 m- A plethora of new phenomena related 
to unique behavior of electrons in graphene were pre¬ 
dicted [HUS]. Although substantial progress has been 
made in getting closer to graphene’s intrinsic behav¬ 
ior mm, there are certain properties that are hard to 
measure in conventional experiments, and consequently a 
complementary approach of decoupling of graphene from 
the environment by levitating small graphene samples 
was proposed HU. This experimental setup has the ad¬ 
vantage of allowing the study of free graphene membrane 
motion at ultrahigh rotation speeds, which can be used 
to apply strain to the membrane mm- Thermal prop¬ 
erties of graphene near its melting point, predicted to 
be close to temperature T = 4510 K [19], are also out 
of reach for standard measurements while easily achiev¬ 
able by laser heating [6]. Finally, the levitation approach 
can be used to conduct research on ultra high vacuum 
(UHV) crystal growth, graphene manipulation and de¬ 
position onto a substrate as well as flakes’ functionaliza¬ 
tion H3HQ. 

In previous experiments m, a quadrupole ion trap 
was used to levitate graphene nano-platelets separated by 
a distance of about 2 mm from the trap electrodes. Dur¬ 
ing the experiment, the time dependence of the center 
of mass motion of the particle and its visible intensity 
were measured at different pressure ( p ) levels, and cir¬ 
cularly polarized light was used to provide spin to the 
particle. While these measurements produced estimates 
of the mass, thickness, and charge, it was also found that 
the time particles stay trapped before they escape was 
limited to a few hours at p < 10“ 7 Torr m However, 
lower p values are required to obtain a low contamination 
environment necessary for material deposition, thermal 
experiments, and high spinning frequency measurements 
on graphene flakes. Consequently, a solution to the prob¬ 


lem of low trap life times is required. 

A similar problem has been known to exist for 
optically trapped microscopic particles, and paramet¬ 
ric feedback cooling was developed as a possible solu¬ 
tion [22] . Experiments conducted in Ref. [22] involved 
silica nanospheres trapped in an optical trap and cool¬ 
ing to 50 mK was observed. In this letter, we discuss 
our implementation of the parametric feedback approach 
for graphene flakes in the quadrupole ion trap, which re¬ 
sults in improvement of the trapping times of the flakes 
compared to the non-cooled case, as well as flakes’ stabi¬ 
lization and cooling in high vacuum conditions. 

The experimental setup used for our measurements is 
similar to the design used in Ref. nzi with modifications 
related to transfer of the particles after trapping and the 
way their motion is detected. The system consists of two 
chambers separated by a gate valve, with ionized parti¬ 
cles introduced into the loading chamber by use of the 
electrospray technique [23] . After trapping in the load¬ 
ing chamber, particles are transferred to the high vacuum 
chamber (see Fig. |T]b) , where the cooldown experiment is 
conducted. The use of two chambers minimizes contami¬ 
nation of the high vacuum chamber by volatile substances 
that are created during electrospray operation. 

When trapped in the high vacuum chamber, charged 
particles are levitated in an electric quadrupole trap[24] 
that consists of two coaxial conical electrodes (Fig. Ilia). 
The inner electrode of the trap is kept at zero voltage 
while an AC voltage V^ c = 300 V at frequency UI/2tt 
(chosen to allow stable trapping of particles) is applied 
to the outer electrode. The AC voltage is created by a 
high voltage amplifier with a signal generator connected 
to its input. Near the trap minimum, the trap potential 
is approximated by [24] : 
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where x,y,z are distances from the minimum along trap 
axes X , F, F (See Fig. [lja), and a Xj y jZ and zq are param¬ 
eters dependent on the trap geometry. For our system 
zq = 2.19 mm as follows from electrostatic simulations, 
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Figure 1. Diagram of experimental apparatus: (a) The par¬ 
ticle trap consists of an outer electrode with a slot (red) and 
an inner electrode (green); (b) Schematic setup of the experi¬ 
ment as viewed from above the trap. Inside the high vacuum 
chamber auxiliary electrodes (Vi, V2, V3) are used to null stray 
electric field; (c) View of the prism mirror beam splitter as 
seen along the direction of the light scattered from the par¬ 
ticle. X ', Y' and Z' show the projections of the main axis of 
motion of the particle onto the lens’ focal plane for three dif¬ 
ferent orientations of the prism beam splitter. Photodectors 
are shown schematically as A and B. 


which are described in full detail in Ref. [25] . 

It can be seen from symmetry considerations (See 
Fig. [I] a) that motion in the X and Y directions is non¬ 
degenerate only if there is some intentional trap asymme¬ 
try introduced to the system. This asymmetry is crucial 
for our experiments, since we rely on the ability to re¬ 
solve all three degrees of translational motion present in 
the optical signal coming from a single detector. The 
splitting of the degeneracy is done by making a slot in 
the outer electrode, which allows us to distinguish mo¬ 
tion along each of the trap axes by its frequency. 

In the pseudopotential approximation, the particle 
motion has three main frequencies of oscillation m : 
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where q and m are the charge and mass of the parti¬ 
cle. The pseudopotential approach is applicable when 
cjy, uj z « £2, conditions which are true for our ex¬ 
perimental parameters of u) Xi y jZ /2ir and £2/27r which lie 
in the range of 150-1000 Hz and 15-30 kHz respectively. 

Charge to mass ratio values of about 10-100 C • kg -1 
are typical for particles trapped using the electrospray 
technique, and the data presented in the paper are taken 
on one of these particles. Special care is taken to make 
sure that only one particle is transferred and studied in 
the high vacuum chamber. When multiple particles are 
trapped, the parameters of the trap (voltage and fre¬ 


quency) are varied towards the unstable trapping con¬ 
ditions, resulting in the consecutive expulsion of the par¬ 
ticles from the trap. The process is repeated until a single 
particle is left over, the case that is easily distinguishable 
visually via the camera. To minimize charge loss of the 
particle, laser power for all experiments was limited to 1 
mW or lower. Discharge observed at higher laser pow¬ 
ers may be related to evaporation of the solution residue 
left on particles after their preparation [26]. A 671 nm 
circularly polarized laser propagating along the Z axis of 
the trap was used to stabilize the rotation of the parti¬ 
cle, and minimize brightness fluctuations attributable to 
nanoflake reorientation m- In the presence of the cir¬ 
cularly polarized laser, the flake lies in the X — Y plane 
and is spinning around the axis parallel to the Z axis of 
the trap. 

Motion of the trapped particle was tracked using a 
532 nm linearly polarized laser beam (E || Y) with a 
width of 0.16 mm (See Fig. [l]b). The light scattered 
from the particle is focused by a lens with a 40 mm di¬ 
ameter located at a distance of 12.5 cm from the center of 
the chamber and is afterwards divided by a 50:50 beam 
splitter. One of the split beams is used for rough particle 
positioning with a CCD (charge-coupled device) camera 
and provides information about average position of the 
particle and amount of scattered light at a frame rate 
of 1 sec -1 . The second beam falls on a knifeedge of a 
90° prism mirror beam splitter at the lens’ focal point. 
This directs light to one of two avalanche photodiodes (A 
and B in Fig. |T]c) , depending on the position of the par¬ 
ticle. The difference between the photodetector signals, 
S = ^qr§, can be used to study the center of mass motion 
of the trapped particle. Here we intentionally choose to 
use the normalized difference of signals in order to min¬ 
imize the signal dependence on the laser power and on 
variations of intensity of scattered light. 

Only components of motion that have non-zero pro¬ 
jection onto the image plane in a direction perpendicular 
to the prism knifeedge will contribute to the signal S. 
Indeed, orienting the prism edge orthogonally to the Z 
axis of the trap (see Fig. [ljc, where prism is shown as 
viewed along the deflected light direction from the view¬ 
point at the center of the trap) will maximize the amount 
of motion detected from the particle oscillations in the Z 
direction. Similarly, rotating it by 90 degrees to the di¬ 
rection parallel to the Z axis will zero the signal from 
the oscillations along Z and maximize the X and Y pro¬ 
jections, albeit at amplitude l/y/2 lower than amplitude 
along Z in the previous case. A rotation angle of ^54.7 
degrees is chosen in order to make scaling factors for con¬ 
tributions from each of the frequencies to be the same. 
In this configuration, the fast Fourier transform (FFT) 
of the signal S contains all three peaks, one for each of 
the eigenmodes of motion (Figj2|. 

The signal S is sent to an FPGA (Field- 
Programmable Gate Array) for digital filtering and pro¬ 
cessing as well as for generation of slow feedback (cut-off 
frequency of 4.4 Hz) that is used to keep A — B by ad- 
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Figure 2. Power spectrum of signal S at (a) p = 4 mTorr 
and (b) p — 4 x 10 -7 Torr. 


justment of lens’ position on a piezostage. The speed 
of this feedback is set to be much slower than the typi¬ 
cal eigenfrequencies of the particle motion to avoid any 
cross-talk between the feedback and the motion. 

The second and main purpose of the FPGA is to gen¬ 
erate a feedback signal that is fed to the signal generator 
as an amplitude modulation (AM) input. The feedback 
signal Sf ee d = gSS is created by multiplying the filtered 
signal S by its derivative [22], where g is the gain setting 
of the feedback. For sinusoidal signals, this operation cre¬ 
ates an output at twice the particle oscillation frequency. 
The feedback modulation is applied to the trap potential 
Vac = Va C (1 + Gam Sfeed), where Gam is the AM value 
set at the signal generator. 

Once the p is lowered below 1 mTorr, the cooldown 
of the particle motion is observed if proper feedback con¬ 
ditions are chosen. We have used AM setting of 30% at 
feedback gain of g = 2.4 x 10 -4 s to show the cooling de¬ 
pendence on p. At these parameters, all three degrees of 
particle motion (FigjSJa) reach T^20Katp = 4x 10 -7 
Torr. In the region between 10“ 3 and 10“ 7 Torr, T is 
strongly affected by the pressure, and even changing p 
from 10 -6 Torr down to 4 x 10“ 7 Torr improves cooling 
from T ~ 40 K down to T = 20 K. 

To prove that the AM setting of 30% chosen for the 
experiment is optimal and to show that parametric feed¬ 
back leads to cooling, we measure the T dependence on 
the AM setting. It can be seen in Fig|3]b that higher 
AM values result in better cooling and there is a tran¬ 
sient region, where cooling efficiency is different for each 
of the degrees of freedom. It is not currently clear why 
this is the case and more experiments are required before 
making any statements about the source of this behav¬ 
ior. At the same time, one can notice that increasing the 
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Figure 3. Dependence of T of the center of mass motion 
along the trap axes X (red), Y (green) and Z (blue) on: (a) 
p ; (b) amplitude modulation Gam [at p = 4 x 10 -7 Torr]; (c) 
voltage V 3 (values plotted relative to the nulling setting). 


AM setting to values higher than 30% does not improve 
the cooling, which makes Gam = 30% the best value 
for experiments. Absence of further cooling for larger 
AM values indicates that noise present in the feedback is 
likely limiting cooling. 

The sensitivit*y of cooling to the feedback noise is 
enhanced by any non-zero DC electric fields that can be 
present in the chamber (arising, for example, from patch 
potentials PQ EZD- This makes it necessary to balance 
out their contribution by nulling the total DC field act¬ 
ing on the particle. The cooling optimization is done by 
finding the settings on auxiliary electrodes (see Fig. 0b) 
that correspond to minimum T of the particle. Although 
these settings tend to change with time the scale of 
their drift is on the order of days. Therefore, their val¬ 
ues need to be found only once during the measurements. 
Typical dependence of cooling on electrodes’ voltage in 
the vicinity of the nulling settings is presented in Figj3]c. 
For simplicity of discussion, the figure shows what hap¬ 
pens when only voltage Vs is offset from its nulling value, 
but similar voltage dependencies are observed for the two 
other electrodes. The presence of a distinct minimum in 
the T data proves that we have indeed found the volt¬ 
age settings that minimize the effect of DC fields on the 
particle motion. The cooling Ts at the minimum are 
T x = 14 ± 8 K, T y = 7 ± 4 K and T z = 16 =b 5 K for X, 
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Figure 4. Time dependence of the particle motion frequency 
for Z direction. Quantized charge jumps found by fitting the 
jump step to give integer jump values are shown by dashed 
lines [N corresponding to each jump is shown on the right]. 


Y and Z degrees of motion correspondingly. 

Additional characterization of the flake, including di¬ 
rect measurement of T, can be done if we can extract 
its mass and charge values. The estimation of these pa¬ 
rameters is complicated since they usually enter physical 
equations as a ratio (e.g. see Eq. ([ 2 ])). The coupling of 
charge and mass can be overcome if some process that 
affects only charge or only mass exists (or is created) in 
the system. One of the measurements allowing this is the 
discharge of the particle that is observed if large enough 
laser power is applied. 

The discharge sequence at laser power of 1.2 mW 
with a distinct discharge quantization is shown on Fig. [4] 
While only dependence of frequency of motion along the 
Z axis on time is presented, jumps in frequency consistent 
with the depicted ones are observed for motion along X 
and Y axes. We can find Slj that corresponds to the min¬ 
imal possible jump (when charge of the particle changes 
from N e to (N =b 1) e) by fitting each of the frequency 
jumps to be equal to N5w, where N is integer. As follows 
from Eq. mass of the particle can be expressed via 
5lj: 


^ 


n 


AC 


(3) 


where Sq = e j28H3Q] . This expression gives us an esti¬ 
mate of m ~ 1.4 x Hr 17 kg for the particle we present 
data on. Total charge q « +970 e is found from the ratio 
q/Sq = oj z /S(jJ z , where u^/27r = 730 Hz at the beginning 
of discharge. It should be pointed that the estimates are 


given assuming that only charge is changing during the 
frequency jumps. This assumption is based on the fre¬ 
quency behavior expected from Eq. ([2|. In high vacuum 
we expect mass to get only smaller since possible changes 
are limited to desorption of volatiles stuck to the particle 
surface and to falling off of loosely attached layers. The 
resulting increase of the eigenfrequency would be con¬ 
trary to the observed time dependence. 

Further improvement in the cooling and control of 
the particles can be gained by using phase-locked loops 
(PLL) for separate cooling of each degree of motion. This 
narrows the bandwidth of the noise that can couple to 
the feedback during the feedback generation. While the 
bandwidth of the current experiment is limited by the 
necessity to keep the signal from all eigenfrequencies, for 
PLL the bandwidth of the loop can be made as small as 
the width of an individual eigenfrequency peak [31]. In 
high vacuum conditions this leads to bandwidth smaller 
than 1 Hz and a drastic improvement of the signal to 
noise ratio is predicted if PLL feedback is implemented. 
PLL will also eliminate noise coupling related to the 
presence of nonzero electric fields. However, it should 
be pointed out that the PLL feedback will have its own 
challenges, since every time a discharge event occurs, the 
frequency of the oscillations may move out of the PLL 
bandwidth range, and during the required reset of PLL 
a temporary lack of proper cooling can lead to the loss 
of the particle. 

The ability to extract the mass of the particle through 
its discharge allows for a separate estimate of particle’s 
T if a proper calibration of signal S in terms of spatial 
displacement is done. However, it should be pointed out 
that there is an intrinsic nonlinearity in our signal detec¬ 
tion that makes the T extraction potentially inaccurate 
and proper calibration is dependent on exact value of the 
reflectivity of the particle as well as high signal to noise 
ratio. 

In conclusion, we have presented experimental results 
on cooling and stabilization of the graphene nanoplatelets 
in high vacuum conditions. It was shown that at p of 
4 x 10 -7 Torr the particle can stay trapped on a time 
scale longer than a few days and cooled down below 20 
K in all three degrees of motion, reaching particle local¬ 
ization below 1 pm in the cooled state. Achieving stable 
motion of the particle allows us to observe the discharge 
of the particle and to extract its mass and charge. Fi¬ 
nally, it was shown that the elimination of stray fields 
was necessary to achieve optimal cooling. 
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